Magnetic Properties of the Metamagnet Ising Model in a three-dimensional Lattice in 

a Random and Uniform Field 
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By employing the Monte Carlo technique we study the behavior of Metamagnet Ising Model in a 
random field. The phase diagram is obtained by using the algorithm of Glaubr in a cubic lattice of 
linear size L with values ranging from 16 to 42 and with periodic boundary conditions. 
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I. INTRODUCTION 

Random fields and disordered magnetic systems have 
been a considerable source of research in recent years. [l| 
The random field Ising model (RFIM) has been one of 
the most interesting subject of research in Physics of 
Condensed Matter in the last fifteen years and it oc- 
cupies prominence position when we deal with disorde- 
red systems. 0, Q The types most common of disorder 
are represented by i) disorder in the bonds and ii) ran- 
domness in the strength of the applied magnetic field. 
In this model the disorder depends on the applied ex- 
ternal magnetic field and although RFIM has deserved 
many investigations from both experimental and theo- 
retical points of vieWjQl no conclusive result has been 
achieved for the understanding the nature of the phase 
transitions and critical behavior. On the other hand, 
questions as the lower critical dimension 0] and the 
existence of a static phase transition have already been 
solved from the theoretical point of view. However, ques- 
tions as the existence of the tricritical point arc still 
opened. [7] The relevance of RFIM is due to the fact that 
it is the simplest to describe the essential physics of vari- 
ous class of disordered systems, which includes: (i) struc- 
tural phase transitions in random alloys, @ (ii) com- 
mensurate charge-density-wave systems with impurity 
pinning, llCl (iii) binary fluid mixtures in random po- 
rous media, [121 (iv) mel ting of intercalates in layered 
compounds such as TiS2,[i3 (v) frustration introduced 
by the disorder in interacting many body systems, be- 
sides explaining several asp ects of electronic transport 
in disordered insulators |14| and (vi) systems near the 
metal-insulator transition. jl5l [T6|ln the last years, the 
physics of the hysteresis, of the avalanche behavior and 
of the origin of self-organized criticality[l^ have been 
modelled employing the non-equilibrium behavior of the 
RFIM. In particular, a new class of problems, such as 
self-generated glassy behabior, has been studied through 
the non-disordered model with infinitesimal random fi- 
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eld. [l7| Recently, random magnetic fields have been ap- 
plied also in metamagnet systems like Ising model and 
importants results have been obtained. [18, 19] An ideal 
metamagnet crystal can be constructed overlapping iden- 
tical layers of spin with ferromagnetic coupling between 
spins next neighbors of each layer and with antiferro- 
magnetic coupling next neighbors spins of the adjacent 
layers. One another way to think a metamagnet crystal 
is to consider a cubic crystal with the couplings between 
the first neighbors as being of antiferromagnetic charac- 
ter and the couplings between the second neighbors as 
being ferromagnetic. In these systems, only the compe- 
tition between the ferromagnetic and antiferromagnetic 
ordering are interesting. However, the application of a 
random and uniform magnetic field can yield the appe- 
arance of new phenomena and a richer critical behabior 
becomes possible and, in particular, in the present work 
we study the RFIM applied to a metamagnet system in 
a cubic lattice by employing Monte Carlo Method. 



II. MODEL AND SIMULATION 

Metamagnet ideal can be considered as a set of 
spins with uniaxial anisotropy ferromagnetic interactions 
within each (Jl > 0) and antiferromagnetic interaction 
(J2 < 0). The hamiltonian model for spin-1/2 is given 
by: 

N 

"H = - ^ JiaiOj - ^ J2C^i(Jk - '^{h - hi)(Ji , (1) 



where the first is executed on all pairs of spin nearest- 
neighbors on same plane and second sum run over all 
pairs of spin nearest-neighbor in parallel planes, h is the 
strength of the external uniform magnetic field and hi 
is the random magnetic field which obeys the bimodal 
distribution given by: 
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P{h,) = -[6ih, - hr) + 6ih, + hr)] 



(2) 



where hr is the strength of the random field. 

To study this system we employed Monte Carlo simu- 
lation technique [20| by using the algorithm of Glauber 
in a cubic lattice of linear size L with values ranging from 
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16 to 42 and with periodic boundary conditions. To re- 
ach the equilibrium state we take, for guarantee, at least 
2 X 10** Monte Carlo steps (MCs) for all the lattice si- 
tes we studied and more 3 x 10* MCS to estimate the 
average values of the quantities of interest. In our work 
we consider one MCs equivalent trials for change the 
state of a spin of the lattice. 

We calculated the sublattice magnetization per spin 
belong to the different planes by using 
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(3) 
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The transition lines of the phase diagram were obtai- 
ned from the staggered magnetization, and magne- 
tization m. Thus we calculate — uia — ttib and 
771 = ruA + 1^13 that are our parameters of order for anti- 
ferromagnetic and ferromagnetic phase, respectively. In 
the above equations [ • • • ] denotes the average over the 
disorder and ( ■ ■ • ) denotes the thermal average. 



III. RESULTS AND CONCLUSIONS 

The Figure 1 illustrate the complete phase diagram, for 
two selected values of the random field hr, showing the 
continuous and discontinuous transition lines separating 
the antiferromagnetic and paramagnetic phases. The tri- 
critical point, which is indicated by an open square, joins 
these two lines. The points of the phase diagram are ob- 
tainded from the knowledge of the point of maximum of 
the curve for the susceptibility, see Figure 2. It was ob- 
tained of the results by fixing H and changing T. In the 
case of the discontinuous transition, we determined the 
magnetization curve as a function of the field, for a fixed 
value of temperature Figure . This procedure is not an 
efficient to localize the tricritical point one, because it is 
difficult to distinguish a continuous from a discontinuous 
curve, especially near the tricritical point, but it gives an 
idea of the range of values of the field where the tran- 
sition is of first order. In this work the location of the 
tricritical point was achieved through the disappearance 
of hysteresis [HI . For a fixed value of temperature, we 
drew the magnetization curves for increasing and decre- 
asing values of the magnetic field. 

(In Figs., we show these curves for a system of size 
L = 40, and for three values of temperature near the 
tricritical point. Our estimate for the tricritical tempe- 
rature is T=.) 

On the other hand, Figure 3 shows the behavior of the 
magnetization m and the staggered magnetization rris as 
function of the temperature for different values of the 
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Figura 1: Phase diagram of metamagnet in a cubic 
lattice in the plane t — h for hr — and hr — I. 
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Figura 2: Susceptibility curve versus temperature, 
KbT/ J, for difi'erent values of external field. 



field external uniform. In this diagram we can observe as 
the random field aff'ects the behavior of the system. 

The behavior of the magnetization when the exter- 
nal field varies for different values of the temperature is 
shown in Figure 4. For T = 2 one get a first-order tran- 
sistion and for T = 3 and T = 4 a transistion of second 
order. 

We can observe in Figure 5 the variation of the suscep- 
tibility with the external field for different temperatures. 
In Figure 6 we can see better as the random field affects 
the order of the system, in this Figure we kept the tem- 
perature and we calculate the magnetization with the va- 
riation of the field uniform for different values of random 
field. 

In summary, the present Monte Carlo simulations for 
a metamegnet Ising model in a random and uniform field 
show that the phase diagram in the plane uniforme field 
h versus temperature present continuous and first-order 
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transition lines separated by possible tricritical points. 
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Figura 3: Curves of magnetization m and staggered 
magnetization rUg for different values of external field. 
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Figura 5: Curves of magnetization m versus external 
field for different values of temperature. 




Figura 6: Curves of magnetization m versus external 
field for different values of random field hr- 



